FINITE DIFFERENCE SCHEMES FOR THE SYMMETRIC 
KEYFITZ-KRANZER SYSTEM 
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Abstract. We are concerned with the convergence of numerical schemes for 
the initial value problem associated to the Keyfitz-Kranzer system of equations. 
This system is a toy model for several important models such as in elasticity 
theory, magnetohydrodynamics, and enhanced oil recovery. In this paper we 
prove the convergence of three difference schemes. Two of these schemes is 
shown to converge to the unique entropy solution. Finally, the convergence is 
illustratred by several examples. 
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1. Introduction 

In this paper, we consider difference methods for the Cauchy problem for the 
n X n symmetric system of Keyfitz-Kranzer type 



(1.1) 



]u{x,0)=uo{x), X S M, 

where T > is fixed, u = (u^^), . . : E x [0, r) M" is the unknown vector 
map with |u| — \/ w^^)^ + • • • + — ^Uq"'^'', . . . , Wq"''^ the initial data, and 

(/) : M — > E is given (suflEicicntly smooth) scalar function (see Section 2 for the 
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complete list of assumptions). Systems of this type was first considered in [10, 12] 
and later on by several other authors [ , as a prototypical example of a non-strictly 
hyperbolic system. This type of system is a model system for some phenomena in 
magnetohydrodynamics, elasticity theory and enhanced oil-recovery. This system 
also has similarities to a model of chromatography [j ] and to a model describing 
polymer flooding in porous media [ '-]. Note that since (/) is a function of we 
call (1.1) a symmetric Keyfitz-Kranzer system. A non-symmetric version of the 
Keyfitz-Kranzer system reads 



(1.2) 



Ut + {U(I){U,W1,W2, ■ ■ - jWn))^ = 0, 

{uWi)t + {uWi(l){u,Wi,W2, ■ ■ ■ ,Wn))x = 0, i = 1,2, ... ,71. 



Existence of global bounded weak solutions to (1.2) has been studied by Lu [ ' '] for 
a specific choice of 0. 

For the flux function F{u) — u(f>{\u\), a straightforward calculation shows B{u) — 
dF{u) is the matrix with entries 

Bi,j{u) = <j){\u\)Sij +(j)'{\u\)^-^, i,j ^1,2, ■■■ ,n, 
where dij is the Kronecker delta, given by 




The matrix B{u) is symmetric, therefore the system (1.1) is hyperbolic, that is, all 
the eigenvalues of B{u) are real and the corresponding collection of eigenvectors is 
complete. It is easy to see that the first eigenvalue of B{u) is Ai = 0(|u|)-|-(/)'(|m|) |m| 
and other n — 1 eigenvalues are = 0(|u|), i — 1,2, •• • ,n — 1. The presence of 
repeated eigenvalues shows that the system (1.1) is not strictly hyperbolic. 

Due to the nonlinearity, discontinuities in the solution may appear independently 
of the smoothness of the initial data and weak solution must be sought. A weak 
solution is defined as follows: 

Definition 1.1. We say u{x,t) a weak solution to (1.1) if 

D.l u{x,t) £ L°°(IR X M+). 

D.2 For all test functions G C(f (M x [0,oo)) 

(1.3) // utpt + u(j){\u\)ipx dxdt + / UQtp{x,0) dx = 0, 
J Jmxm+ Jr 

It is well known that weak solutions may be discontinuous and they are not 
uniquely determined by their initial data. Consequently, an entropy condition must 
be imposed to single out the physically correct solution. Therefore the Cauchy 
problem is viewed in the framework of entropy solutions. For (1.1), an entropy 
formulation was first introduced by Freistiihler [5, (i], and independently, by Panov 
[!()]. An entropy solution to (1.1) is defined as follows: 

Definition 1.2. A bounded measurable function u{x, t) is called an entropy solution 
to (1.1) if 

D.l For all test functions e C^f (M x [0, oo)) 

(1.4) // uipt + u(j){\u\)ipx dxdt + / UQtp{x,0) dx = 0, 
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D.2 r — \u\ is an entropy solution (in the sense of Kruzkov [ ]) of the scalar 
conservation law 



(1.5) 




uo{x) \ . 



Regarding the existence, uniqueness of solutions and continuous dependence of 
solutions on the initial data we have the following result 

Theorem 1.1. The system (1.1) has the following properties: 
(E) The system has a solution for any uq G L°°{M.). 

(U) For such uq, there is precisely one solution u with the property that r — \u\ 
satisfies the scalar conservation laws (1.5) and Kruzkov's entropy criterion. 
(S) This solution u depends L\^^{M.) continuously on the initial data uq. 

This theorem was first proved in [(>] by using the famous equivalence result of 
Wagner [19]. The key idea behind this proof is to view the system (1.1) as an 
extended system, consisting of (1.1) and an additional conservation law satisfied by 
r (1.5), with Wagner's transformation theory. On the other hand, in [IG], Panov 
gave a "direct" proof of both existence and uniqueness. The existence was proved 
by showing the convergence of the singularly perturbed problems 

ul + iu'(t){\u^\))^^eul^, 

to an entropy solution as e — > 0. The idea behind the existence proof was first to 
show the existence of a measure-valued solution vit^x) of the Cauchy problem (1.1). 
Then he showed that indeed i'(t,x) is regular: i^jt ^.^(u) = S{u — u{t,x)),u(t,x) S 
L°°{M. X K+,K") and consequently this gives existence of a solution to (1.1). 

In view of the analytic properties of the solutions of (1.1), several different meth- 
ods for computing the solution suggest themselves. Foremost among these methods 
is Glimm's scheme [n]. Regarding other numerical methods, it is tempting to use the 
equation satisfied by r, and view r as an independent variable. Defining v G S*""^ 
hy vr = u, we formally have that 

(1.6) rt + ir^{r))x=0 

(1.7) irv)t + {r(j)ir)v)x = 



or 



(1.8) vt + (l){r)vx ^ 0. 

As a strategy, one can then solve (1.6) first, and then either (1.7) or (1.8). These 
should then hold subject to the constraint \v\ = 1. Without this constraint, (1.6)- 
(1.7) is a "triangular" system of conservation laws, see [.!]. Using any monotone 
scheme for (1.6) and (1.7) will ensure the strong convergence of the approximate 
solutions to (1.6) and the weak-star convergence of the approximate solutions to 
(1.7). This approach was used in [7]. To show that u ^ rv is an entropy solution to 
(1.1), one must show (for the approximations) that |w| = 1 in the limit if \vo\ = 1. 
In this paper, for the semi-discrete scheme, we discretize (1.1) in space and show 
the convergence of approximate solution to a weak solution of (1.1). But we are 
unable to extend our analysis to the fully discrete scheme based on discretizing 
(1.1). To overcome this difficulty, we propose another scheme based on discretizing 
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(1.6)-(1.8) and prove the convergence of approximate solution to unique entropy 
solution of (1.1). 

The present paper can be divided into four parts: 

(1) In Section 2, we present the mathematical framework used in this paper. 
In particular, we used a compensated compactness result in the spirit of 
Tartar [1.7] but the proof is based on div-curl lemma and does not rely on 
the Young measure. 

(2) In section 3, we propose an upwind semi-discrete finite difference scheme 
and prove the convergence of the approximate solution to the weak solution 
of (1.1). The idea behind this proof is to prove first the strong convergence 
of approximate solution r\x = |ua2:| using the compensated compactness 
technique [17, 2]. Then prove a B.V. estimate for r = u/{u, e), where e is a 
unit vector in M". Then Kolmogorov's theorem, combined with the strong 
convergence of ta^i gives the strong convergence of approximate solution 

UAx- 

(3) In section 4, for a fully discrete scheme, we are only able to conclude that 
u is only a distributional solution of 

ut + {u(j>{r))x = 0, 

for some r such that |m| < r. We propose another fully discrete scheme 
relying on explicit decoupling of the variables r and v expressed by the 
"nonconservative" formulation (1.6)-(1.8) 

\vt + (t)ir)vx =0, 

with r(0) = 1^(0)1. It is not difficult to show the convergence of r^x to r, 
r being the unique entropy solution of (1.6), and the strong convergence 
of fAx- In order to conclude that u = rv is the unique entropy solution of 
(1.1), one has to show \v{x, t)| = 1 and this has been achieved in this paper 
using Wagner transformation [l!)] (see Section 2 for more details). 

(4) Finally, in Section 5, we test our numerical schemes and provide some 
numerical results. 

2. Mathematical Framework 

In this section we present some mathematical tools that we shall use in the 
analysis. To start with the basic assumptions on the initial data and the funtion 
(/)(r), we assume that </> is a twice differentiable function (f> : [0, oo) — > [0, oo) so that 

A.l 0(0) = 0, 4'{r) > and (/)'(r) > for ah relevant r; 

A. 2 4>{r), (j)'{r) and (/'"(r) are bounded for all relevant r; 

A.3 meas |r 2(p'(r) + r0"(r) = o| = 0; 

A.4 |uo| e L^{U) n L°^(R) and |uo| G B{K) for any constant K in M, where 
B{K):={f I inf / > Ca-}, and 

xeAiJ.K) 

( xe {-oo,K] I 3e > I 

I with f{y) > liminf /(z), for a.c y G {x — e,x) \ ' 

K z^x J 

Here Ck is a positive constant depending on K; 



A{f,K) :-- 
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A. 5 Mo G r^-, where Ts is the cone 

Ts {u e M" I S\u\< (e,u)} 

for some fixed unit vector (which we without loss of generahty choose as 
e — {1, . . . , l)/^/n), and 6 is a, fixed number in the interval {^y{r^^^T)Jn, 1). 

Next, we recapitulate the results we shall use from the compensated compactness 
method due to Murat and Tartar [J ", 17]. For a nice overview of applications of 
the compensated compactness method to hyperbolic conservation laws, we refer to 
Chen [2] . Let A4 (M) denote the space of bounded Radon measures on E and 



Co(M) = Uje C{R) lim ^P{x) = 

[ |a;|— )-oo 

li fie M{R), then 

{^,iP)= ^Pdn, for aU 7/;eCo(M). 



Recall that fi e M{R) if and only if <C IIV-ILocfR) for aU ^ e Co{R). We 

define 

IImIIm(r) -sup{|(/i,V)| -.^ eCo{M.),U\\L^(R) < !}• 

The space is a Banach space and it is isometrically isomorphic 

to the dual space of ^Co(M), ||-|lioc(K)^ , while we define the space of probablity 
measures 

Prob(IR) = e A^(]R) : ^is nonnegative and |l/i|l_A4(R) = !}• 

Before we state the compensated compactness theorem, we recall the celebrated 
div-curl lemma. 

Lemma 2.1 (div-curl lemma). Let ft be a bounded open subset o/M^. With e > 
denoting a parameter taking its value in a sequence which tends to zero, suppose 

D'' ^ D in{L^{VL)f, E m {L^{n)f, 

{divD^I^^Q lies in a compact subset of H^^{n), 

{curlE'^I^^Q lies in a compact subset o/iJj^^(51). 

Then along a subsequence 

■ E'' D ■ E in V'{n). 

We shall use the following compensated compactness result. 

Theorem 2.1. Let ft C M. x IR+ be a bounded open set, and assume that {w'^} 
is a sequence of uniformly bounded functions such that < M for all e. Also 
assume that f : [~M, AI] — > M is a twice differentiable function. Let u'^ ^ u and 
f{u'^) V, and set 

{r^,{s),q,{s)) = {s-k,f{s)- f{k)), 

(2.1) 



iV2{s),q2{s))=^ns)-fik),J^ {f'{6)Yde 

where k is an arbitrary constant. If 

'ni{'^'^)t + <li{u^)x is in a compact set of H^^^{Q) for i — 1, 2 



6 



U. KOLEY AND N. H. RISEBRO 



then 

(1) V = f{u), a.e. {x,t), 

(2) u" u, a.e. {x, t) i/meas{u | f"{u) = 0} = 0. 

For a proof of this theorem, see the monograph of Lu [ l'']. A feature of the com- 
pensated compactness result above is that it avoids the use of the Young measure 
by following an approach developed by Chen and Lu [13, 2] for the standard scalar 
conservation law. This is preferable as the fundamental theorem of Young measures 
applies most easily to functions that are continuous in all variables. 

The following compactness interpolation result (known as Murat's lemma [15]) 
is useful in obtaining the H^^^ compactness needed in Theorem 2.1. 

Lemma 2.2. Let Q. he a hounded open subset ofM.'^. Suppose that the sequence 
{£e}^^Q of distributions is bounded in W~^''^{fl). Suppose also that 

where is in a compact subset of H~^{fl) and {£2,e}£>o ^ bounded 

subset of Ai\oc{^)- Then {/^gj^^g is in a compact subset o/i/[~^(Sl). 

Next, we shall need Kolmogorov's compactness lemma. 

Lemma 2.3 {Lj^^ compactness, see [ ]). Let u'^ : M. x [0, oo) — > M &e a family of 
functions such that for each positive T , 

\u'{x,t)\ < Ct, {x,t) e M X [o,r] 

for a constant Ct independent of e. Assume in addition that for all compact i? C M 
and for t G [0, T] 



sup / \u''{x + ^,t) - u''{x,t)\ dx <vb.t{\p\)-, 
\e.\<\p\JB 

for a modulus of continuity vb,t- Furthermore, assume for s and t in [0,T] that 
\u'^{x, t) — u^{x, s)\ dx < u!B.Ti\t — s\) as e 10, 



for some modulus of continuity ujb.t- Then there exists a sequence {sj} such that 
for each t G [0,T] the function {u'^^(t)} converges to a function u{t) in L^^^iM.). 
The convergence is in C{[0,T]; Ll^^{M.)). 

Finally, we state the following result related to Wagnar transformation theory. 

Lemma 2.4 (Wagner Transformation, see [19, 5]). For any n G N, there is a one-to- 
one correspondence between (equivalence classes of) bounded Lebesgue measurable 
solutions (r,rv) : — > [0, oo) x M™ to the system (1.6)-(1.7) which satisfy 

rO poo 

r{x,t)dx^ / r{x,t)dx ~oo 



and (equivalence classes of) weak solutions {t,v) to the system 
(2.2) 



in which t is a Radon measure in M'^ which dominates Lebesgue (outer) measure 
A2 (i.e., T > k\2 for some k > 0), f is the density of the absolutely continuous part 



KEYFITZ-KRANZER 



7 



of T with respect to X2 , and w : ^ 

hounded and Lehesgue measurable. 
This correspondence is established through transformations T : {x,t) — )■ (y{x,t),t) 
defined relative to any hounded measurable solutions to (l.C) by 

(2.3) ^{x,t)^r{x,t), ^{x,t) = ~cl,{r{x,t))r{x,t), y(0,0) = 0, 
ox at 

namely setting 

(2-4) , ' ' 

V — V o 1 

Observe that using (2.2), (2.3) and (2.4), it is easy to conclude that \v\ = 1. 



3. A SEMI-DISCRETE FINITE DIFFERENCE SCHEME 

We start by introducing some notation needed to define the semi-discrete finite 
difference schemes. Throughout this paper we reserve Ax to denote a smah pos- 
itive number that represent the spatial discretization parameter of the numerical 
schemes. Given Ax > 0, we set Xj — j Ax for j £ 'Z and for any function u = u(x) 
admitting pointvalues we write Uj = u(xj). Furthermore, let us introduce the 
spatial grid cells 

where Xj±i/2 — Xj ± Ax/2. Let D± denote the discrete forward and backward 
differences, i.e.. 

The discrete Leibnitz rule is given by 

D±{ujVj) — UjD±Vj + Vj±iD±Uj 

Furthermore, for any function /, using the Taylor expansion on the sequence 
f{uj) we obtain 

D±fiu,) = nu,)D±u, ± ^f"{i,±^){D±u,)\ 

for some between Uj±i and Uj. We will make frequent use of this, which 

states that a discrete chain rule holds up to an error term of order Ax{D±Uj)'^ . To 
a sequence {uj}j^^ we associate the function defined by 

UAx{x) = y^Mjlj, (x). 

Similarly, we also define rAa; as 

rAx{x) = ^rjl/^(x), 

where 1a denotes the characteristic function of the set A. Throughout this paper 
we use the notations u^.^, r^x to denote the functions associated with the sequence 
{■f^jljgz s-^^d {''jljez respectively. For later use, recall that the L°°(IR) norm, the 
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i^(M) norm, the L'^(M) norm, and the BV{M.) semi-norm of a lattice function uax 
are defined respectively as 



\UAx\\li 



Observe that all the eigenvalues of the system (1.1) are positive by our assumptions. 
Therefore we consider the following semi-discrete upwind finite difference scheme 

(3.1) u'^{t) 4- D_ {(j){rj{t))uj{t)) = 0, for j e Z, i > 0, 
with initial values 

(3.2) uj{0) = ^ / ua{x)dx, 

JXj_i/2 

where rj{t) — \uj{t)\. We have that {uj{t)} .^^^ satisfy the (infinite) system of ordi- 
nary differential equations and it is natural to view (3.1) as an ordinary differential 
equation in L^(IR)", since the piecewise constant structure of uax is preserved by 
the evolution equation (3.1). To show the local (in time) existence and uniqueness 
of differentiable solutions we must show that the right hand side of (3.1) is Lipschitz 
continuous in L^(M)". Set 

F{uAx)j = £>- {(t>{rj)uj) . 
The infinite system of differential equations (3.1) can then be written 

at 

We view F{uax)\x as an element in L'^(IR)". To establish that this system has a 
unique solution (at least locally in time) we show that 

(3.3) \\F{uAx)^^x - F{vAx)^^x\\L2(;Si)^ <1\\uAx-VAx\\l2{K)^ 

for some locally bounded 7 = 7('UAa;, wak) and for a fixed Ax > 0. Set fj — \vj\, 
note that 

V] - I < , - \Uj ~Vj\< \Uj -Vj\. 

Then 

2 / 

\\F{uAx)- F{vAx)\\L2{tLy^ < ^(sup|wj| ||0'||^„ I^A^. - ^=A^|li2(R) 

<j\\uAx-VAx\\l2(R)^ , 
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where we have used Assumption A. 2. Therefore F is locally Lipschitz continuous, 
and there is a r > so that the initial value problem (3.1) has a unique difFerentiable 
solution for t G [0,t), if r < oo, then 

lim||uAx(i)llL2(K)„ = 00. 

We shall proceed to show that the norm remains bounded if it is bounded 
initially, so the solution can be defined up to any time. 

Lemma 3.1. Assume that A.l, A. 2 and A. 4 hold, and let {uj{t)} be defined by 
(3.1), and let rj{t) = \uj(t)\. Then 

||»'Ai;(i)|lLi(R) < ||?'Ai(0)|l^i(R) , 

lkA:.(i)|li2(R) < ||rAx(0)||i2(R) , 

lkAr.(t)|lL<»(R) < ||rAx(0)||ioo(R) . 

Furthermore, there is a constant C , independent of Ax and T , such that 

/ I V (r|-s2)0'(s) ds + Ax V0j_i Ax |L>_Ujn dt< C. 

Proof. Let r; be a differentiable function 77 : M" M, take the inner product of 
(3.1) with Vu'ni'^j) to get 

(3.5) j^ii{u,) + D^{cj>,7^{u,)) 

Ax 

+ [{VuV{uj),Uj) - rj{uj)] D^(j)j + <i^?7j-i/2 {D^Uj,D_Uj) = 0. 

Here (pj = <t>{rj), and d^ry denotes the Hessian matrix of 77, so that 

dlvj~i/2 = dlr]{uj^i/2) 

for some Wj-1/2 between Uj and Wj-i. By a limiting argument, the function ri{u) = 
\u\ can be used. If one approximates by convex smooth functions, this means that 

d 

di' 

Multiplying by Ax and summing over j we get 

(3-7) ||rAa;(i)llLi(R) < II"o|Ili(b) ■ 

Furthermore, if j is such that rj{t) > rj_i(t), then since (j) is non-decreasing we 
get rj{t)(j){rj{t)) > rj_i(t)0(rj_i(t)) i.e., D_{rj(j){rj)) > 0. Hence, from (3.6), 
we see that drj{t)/dt < 0. This shows that < rj{t) < sup^ |uj(0)|. Hence, if 
II |uo| ll^ooj-R) < 00, then rAx is bounded independently of t and Ax. 
Choosing 77(1*) — \u\^ in (3.5) we get 

j^r^it) + D_ {r^cj}^) + rp_<j,^ + 0,_iAx \D_u^f = 0. 

We have that 

1 



(3.6) -rj+D^{r,^{r,)<0. 



(r|0,) + rp.cb, - ^ / (s^is) + s^cj,'{s)) ds + ^ I (r| - s') </.'(s) ds 
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where 

(3.8) g{r) = 2 [ {s(l>{s) + s^(l)'{s))ds. 

Jo 

Using this we find that 

(3-9) |kAx(t)|lL2(R) < III"o||Il2(k) , 

since, by the assumption that cf)' > 0, 

{r]~s^) (/)'(s)ds>0. 



rj-i 



Hence ||uAa;(0llL2(R)'» bounded independently of Ax and t. Therefore, there 
exists a differentiable solution UAx{t) to (3.1) for all t > Q. Furthermore, we have 
the bound 

/ (V / ' (r2-s2)0'(s)ds + Aa;V(^j_iAa;|L>_Ujn dt<C, 
Jo y ^ ^ J 

for some constant C which is independent of t and Ax. 

□ 

Now let i5 be a positive constant, and let e be some unit vector in M". Choose 
ri{u) = max {S \u\ — (e, u), 0} . 
and observe that {Vr](u),u) — rj{u) = 0. Furthermore rj is convex, so that 

d 

dt' 

which implies that 

^?7(t^,(t)) <5]?7(w,(0)). 

i i 

We have that 77(11) = if and only if u is in the cone Tg = {u | S\u\ < (e,u)} 
for some unit vector e. If 9 denotes the angle between e and u, then u G if 
cos(^) > S, thus if (5 < 1 this is a cone in M" and this cone is positively invariant for 
(3.1). Observe that there is no loss of generality in choosing the coordinates such 
that e = (1, . . . , l)/^/n. If 1 > S > {n ~ 1) /n, then the invariant cone is in the 
first 2"-tant in M", so that u^'V) > for aU i > if uq G F^. Since = \uj\, it 
follows that u"j\t) = if and only if rj{t) — 0. 

Therefore, if Uq G {u | |u| < i?} n F^, then UAa;(a;,t) is also in this set. This 
enables us to deduce the weak-* convergence of a subsequence (which we do not 
relabel) of {uax}ax>o- 

Let now rji{r) and qi{r) be given by (2.1) for i — I, 2. We then have that 

(3.10) j^m{uj) + D^{q,{r,)) + e,,,=0, 

where 

f{r) ^ r(j}{r), 9i(r) = /(r) - /(fc) and 

e,, = Ax {D^u. f -^(l- -^-^ff^^-^A (^_,^.) . 

r 3-1/2 \ H_i/2 J 
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For any vector u, the matrix u® u is defined as {u ® u)ij = UiUj. We shall now 
find an equation satisfied by 772. First observe that 

Multiplying this with f'{uj) we get 

(3.11) + q'2{r,)D^r, - fir,)^f {r,_,/,) [D^r.f + /'(r,>i,, = 0. 

Set 



^2,. = ^/" b-1/2) [D^r.f. 



Then (3. If) can be rewritten as 

(3.12) j^m{r,) + D^q^r,) + ^q'^{r.,^„^) (D.r.f - f'{r,) (e^,, - e,,,) = 0. 
Finally set 

eaj = ^92(^^-1/2) [D-rj) , 

and 

ei{x,t) = ejj(t) for x G (0:^-1/2, 2:^+1/2] and i = f,2,3. 

Lemma 3.2. Assume that A.f , A. 2 and A. 4 /loW, f/ien we have that e.^ G A^ioc(]Rx 
[0,T)) fori^l, 2, 3. 

Proof. Set = M X [0,r), and let ^ be a test function in Co(il). Note that from 
(3.7) and (3.f0) it follows that 



ei{x, t) dxdt < C, 

where the constant C does not depend on Ax or T. Since ei > 0, this means that 

|(ei,V)l < jj^meidxdt<C\mL^(n)^ 

and thus ei G A^ioc(f^). To show the same for 62 and 63 observe that 

\D^r,\ < \D_u,\. 
Since (p{r) > 0, (3.4) implies that 

fT 

/ Aa; VAa;|i:)_Uj|^ < C, 
Jo ^ 

for some constant C which is independent of Ax and T. We also have that /' and 
/" are locally bounded, and r^^ is bounded, this means that, for i = 2, 3, 



e. 

n 



^(x,t)dxdt<C Ax^Ax{D_rjf dt< I Ax V Aa; |D_itj |^ < C. 
Jo „ Jo 



Thus also 62 and 63 are in A^ioc(fi). □ 
Observe that. Lemma 3.2 implies that also f'{rj){eij — e2j) is in Miod^)- 
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Lemma 3.3. Assume that A.l, A. 2 and A. 4 hold, let uax be generated by the 
scheme (3.1) and set r/^^ = IuakI- Then 

{r]i{rAx)t + qt{rAx)}Ax>o compact in H{^1, 
where rji and qt are given by (2.1) for i = 1,2. 

Proof. Let i — 1 or i — 2, and ip is a. test function in H^^^{fl). we define 
= {ViiT^ Ax) t + qi{r Ax) x,ll^) 

= (^Yl J (^^Vt(rj)i^ix,t) - q,{rj)ip3;{x,t)^ dx^ dt 

^ fX!^/ ^(??»(?'j)V'(a;,i)da;-gj(rj) AxD_-0(a;j+i/2,t)^ j 



\Y1 J (^^Vtirj)ijix,t) + D_qi{rj)ij{xj_i/2,t)^ dxj dt 
^ (^'y*('^j) iKx,t)dx^ dt 

2 

{ip{xj_i/2,t) - i^{x,t)) D_qi{rj) dx 1 c^t 



E 

=: {Ci^i,ip) + (£2,i, V')- 

By (3.10), (3.12) and Lemma 3.2 we know that £,1 e A^ioc(^l). Regarding £i 2 we 
have 



i(>c2,.,V')i= / E 



2:3+1/2 



ipxiy,t)dyD^qi{rj)dx\ dt 



< 



< 



[ (E/ ^ ^j-i/^i [ {^x{y,t))'^ dy\ \D^q^{rj)\ dx\ dt 



lE^^'^'(£ iMx,t))'dx) U\\^^\D^r,\\ dt 

[ {^^x jy\Mx,t)f dx)"\^x^Y.^D^'-^f)"^^ 



dt 



<\\q-\\^^y/Ax(^f f {'ipx{x,t)f dxdt^^'"^ n AxJ^AxiD^rjfdt 



1/2 



< CVAxM^i^n)- 

Therefore the above estimate shows that £2,1 is compact in H~^{^}). By Lemma 2.2, 
we conclude the sequence {r]i{rAx)t + qti^Ax)} ax>o compact in H^^{fl). 



□ 
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Lemma 3.4. Assume that A.l - A. 4 hold, then there is a subsequence of {Ax} 
(not relabeled) and a function r such that rj^^ — )• r a.e. (x,t) € fl. We have that 
r e L°°([0,r];Li(M)). Furthermore, r satisfies 

frt + /(r), <0, xeR, t>0, 
[r = \uo\ , xeR, t = 0, 

in the distributional sense. 

Proof. The strong convergence of rAx follows from the compensated compactness 
theorem, Theorem 2.1 and the compactness of {rii{r/\x)t + Qi{fAx)x}Ax>Q * ~ 
1, 2, i.e., Lemma 3.3. 

To see that r is a distributional subsolution of the conservation law (1.6), multiply 
(3.6) with a non-negative test function ip and integrate over x and t to obtain 

T 

rAx^Jt + f{rAx)iJxdxdt+ I rAx{0,x)i:{0,x)dx 



pT ^ rXj + i/2 nx+Ax 

- / Z^-^('"j)a~ / / {tpxix,t) - ^px{z,t)) dzdxdt. 

Jo J Jxj_i/2 Jx 

The term on the right tends to as Ax — 0, which shows that r is a subsolution. □ 



Let now the vector r, be defined as 



e= ^(1,...,1), 



{uj,ey / 

if Uj 7^ 0. If Uj — set Tj — Tj+i. Observe that this makes sense since uj*'' = 
only if rj — 0. Indeed, we have 

J^uf\t) + uf{t)D^<^, + 0,_i(<)i?_4') = 0. 

If uf{t) > for t < to and ufita) = then duf /dt{to) < 0. If ^'2i(^o) > then 
duf /dt{tQ) > 0, which is a contradiction. Thus if rjg{to) — 0, then rj{t) — for 
all j < jo and t > to. Thus the definition of r-'-' makes sense. 



First note that 



(D-Uj) {uj,e) - Uj {D_Uj,e) 



{uj,e) {uj-i,e) 
Using this, we find 

d {duj /dt — {duj / dt, e)) {uj, e) — {uj — {uj, e)) (duj /dt, e) 
dt^' ^ iu,,ef 

£*_ {uj4>j) (uj, e) — Uj (£)_ [ujcj)) , e) 
(uj,ef 

(D-Uj) (mj, e) ~ Uj {D^uj,e) 
= ^Vj-i ; ^2 



{uj,e) 

(uj-i,e) 
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Set 



Now Tj satisfies 




ifr, >0, 



ifr,=0. 



(3.13) -r,+\,D^T,^Q. 



d_ 
dt' 

Next, we define for any constant J > the set 

M{J; t) := {j I rj{t) < rj_i(t) and j Ax < J}, 

and 

jeA1(J;t) 

Lemma 3.5. Assume that the assumption A. 4 holds. Then we have that Tj{t) > 
STj{0) > SCj > 0. 

Proof. First let j e A4{J;t). In particular this means that rj{t) < rj_i(t). Using 
scheme (3.1), we have 

j^{uj,e) = -{D^{u.i<j)j),e) 

= -D^(f>j{uj,e) - 0j_i(D_Uj,e). 

Therefore, if {uj,e) is a local minima, i.e., if {uj,e) < (uj_i,e) then -^{uj,e) > 0. 
This concludes the proof since SCj < STj{0) < Srj{0) < (ui(0),e) < {uj{t),e) < 
rj{t). Observe that this proof also implies that the set M{J\t) can't increase in 
time. □ 

Lemma 3.6. If the assumptions A.l - A. 5 hold, and if 



(3.14) r«(-,0) 



<C, 1 = 1, 

s.y.(R) 



for some constant C which is independent of Ax, then there is a subsequence of 
{Ax} (not relabeled) and functions t^''^ in C{[0,T]; Ll^^(M.) such that T^l,{-,t) — >■ 
TW(-,i) m Ll^iR) fort = l,...,n. 

Remark 3.1. //r(-,0) is of bounded variation, or if uq is of bounded variation, 
and satisfy A. 4 and A. 5, then (3.14) holds. 

Proof. Note that Xj > 0, and that Xj is bounded in any compact interval. We only 
want to show that Xj is bounded in any compact interval of M, since Kolmogorov's 
compactness theorem gives the convergence in C([0, T]; Lj'Q^(M)). First observe that 



(uj_i,e) 

"i-i^ ^ 

(uj,e) 



~ 6 r^ 



Now if rj > Tj-i, then clearly Xj is bounded. Again if rj < rj-i, thanks to 
Lemma 3.5, we have Xj is bounded. 
Set = D_Tf\ Then Oj satisfies 

(3.15) ^,+x,^,D_e,+e,D_x,=Q. 
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Let r]a{9) be a smooth approximation to \6\ such that 

riaid) > and hm 7y„(0) = Um (^(0)) = |0| . 

Q— ^0 a— >-0 

We multiply (3.15) by r]'^{Oj) to get an equation satisfied by 'i]a{Oj). Observe that 

+ {0jV'M-VamD_X,. 



Hence 



Now let a — to obtain 
(3.16) 



dt 



\e,\ + D^ {x,\e,\)<{). 



If we multiply this with Ax, sum over j and integrate in t, we find that 



(3.17) 



< 



B.V 



< C 



B.V. 



Note that, since T^l.{-,t) has bounded variation and satisfies (3.13), it is L^^^ Lips- 
chitz continuous in t, that is 



(3.18) 



< sup Aj 

LHl-J,J]) jAxe[-j,j] 



for any compact interval [—J, J] of E. Hence, the above estimates (3.17), (3.18) 
and an application of Kolmogorov's compactness criterion (Lemma 2.3) shows that 
r(') = Mm^x^o continuous in t, with values in ^^^^(IR). In other words, the 

convergence is in C([0,T];Lj', 



loc 



□ 



Now we have the strong convergence of r^^ and of tas- This means that also 
UAx converges strongly to some function u in C([0, T]; LJq^(K)) since we have 



(3.19) 



u — t(u, e) = -r^r. 

\t\ 



Theorem 3.1. Assume that A. 1 - A. 5 and (3.14) hold. LetuAx be defined by (3.1) 
- (3.2). Then there exists a function u in L°° ([0,T]; Ll^^{M.)) and a subsequence 
{Axj} of {Ax} such that u Ax u as Axj -> 0. The function u is a weak solution 
to (1.1). 

Proof. We have already established convergence. 

It remains to show that w is a weak solution. To this end, observe that^ 

D- {uAxjHrAxj)) i^{x,t) dxdt = / UAx^(f>{rAxj)D+ip{x,t) dxdt. 

Jo JS. 

As Axj 0, D+ip — > ipx for any e CQ(r2). This means that m is a weak 
solution. □ 

^Here we "extend" the definition of D_ and to arbitrary functions in the obvious manner. 
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4. Fully Discrete Schemes 

In this section, we propose three different fuUy discrete schemes and show two of 
them converge to the unique entropy solution of (1.1). We start by introducing some 
notations needed to define the fuhy discrete finite difference schemes. We reserve 
At to denote a small positive number that represent the temporal discretization 
parameter of the numerical schemes. For n = 0, 1,-- - ,N, where NAt — T, for 
some fixed time horizon T > 0, we set i" = nAt. For any function v{t), admitting 
pointvalues, we let denote the discrete forward difference operator in the time 
direction, i.e., 

+ ^ ^ At 
Furthermore, we introduce the spatial-temporal grid cells 

As before, to a sequence , „ we associate the function uax defined by 

UAx{x,t) ^ ^ w"l/ji(a;,t), 

similarly, we also define r^x as 

jez,Ti>o 

First, we consider the following fully discrete finite difference scheme 
(4.1) Dlu]+D_{u^c^{\u^;\))^0, 
with initial values 



,0 



uq{x) dx. 



^ Ax 

We start by proving the following lemma. 
Lemma 4.1. Suppose uq E L'^{R)" D L°°(M)", and 



X < min ■ 



1 



1/1 



^0 (1 + j|uo||i<»(R)) 



hold. Then 



||"A2;(-,tn)|li2(jjw < ||uo||^2(K)„ , 

(4.2) 2 2 
||uAK(-,ir!.)|lLoo(R)„ < ||uo|Il°°(R)'« ' 

for all n > 0, furthermore 

(4.3) AtAx J2 ^2;|i?_u;f <2||uo|li2(R)„, 



where X = ^ and NAt = T. 
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Proof. For any (differentiable) function rj : M" — >■ M, we can take the inner product 
of (4.1) with V„?7(u") and obtain 

(4.4) DXv{u])+D^{rj{u'l)<P^) 

As before, we choose 'q{u) — \u\^ to get 



The upper hue in the above formula can be rewritten as 

D\ \u-\' + (g-) + ^ [{r^y - s') 0'(.) ds, 

where g is defined in (3.8). In order to balance the two last terms we proceed as 
follows. 



> cl}{0)Ax\^\D^u]\ 

- A (^((1 + \u]_,\) max0'(r))' \D_u]\'' + |uj |' max (/.'(r) 

> (t>{0)Ax (l - C^X (1 + , 

where the constant only depends on 0'. We have the obvious inequality 

n|2 1 2 

rj I - ^ ll"A2:(-,in)|jL2(R)„ • 

We shall not use this, instead we assume that A is so small that 

'3 



(4.5) 1-C^A ( 1 + supr" ) > ^ 



If this holds, setting R = r^, then we have that 



2 



^(.9(v^)-.9(^/^ 



^i-l/2 



^"-A/'(r;Li/2)(i?"-i?^i) 

(l - A/'(r;.Li/2) ^" + A/'(r;.Li/2)^"-i, 
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where /(r) — r(j){r). Hence if 

(4.6) A||/'||^^<1, 

then i?"^^ is dominated by a convex combination of i?" and Rj-i, and thus 
sup,- < sup^- {\u]\}. Therefore, if we assume that 

(4.7) A < min ■ 



we have that r" < r'j < |1mo|Il°°(r) " J- Hereafter, (4.7) is always 

assumed to hold. 

Under the CFL-condition, (4.7), the equation for i?" can be written 
D'^R'J + D^G{RJ) + ^ \D-u']\^ < 0, 
where G(i?) = g{VR). Summing this over n = 0, . . . , iV — 1 yields 

Aa;^i?f + AxAt^^AxlD-i*;'!^ < Ax^i?° < / \uof dx. 

'J n,j 3 "'K 

This finishes the proof of the lemma. □ 

Rewriting yet again the scheme for i?, we have 
(4.8) D^R] + D^G (R^) = el^ + el^, 

where 



= -07_iAa; + Ai |' • 

Let us define 

ei,Ai^ = e"^ for x G (xj„i/2, Xj+i/a] and t e [t„,i„+i), and i = 1,2. 

Lemma 4.2. Assume that A.l, A. 2 anrf A. 4 /loZd, f/ien we have that Ci^Ax G 
Mioc(f^) /orz = 1,2. 

Proof. First we observe that there is a constant C (independent of n, j and Ax 
such that 

\el^ \ < CAx\D_u]\'^ . 

We also have that 

\el, I < £ {r- - r-_,f {r^ + r^,) < CAx \D^u]\' , 
for some constant C depending of (f) and HuoH^oo^u-)- Then, for a test function 



for 1 = 1,2. This finishes the proof. □ 
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Multiplying (4.8) by G'(i?;0 we find 



el 



(4.9) D\G^^ + D^^l^ = G' [K]) [el^ ^ o,^^ 



where 



We find that 



q2{R) = / {G'{e)Y de. 



G"iR)^lrir) + ^. 
2 r 

This is potentially unbounded as r | 0, hence we shall assume that 4> is such that 

(l>'(r) 

(4.10) lini--^<cx). 

riO r 

We are going to show that also the right hand side of (4.9) is in A4\oc{^)- First 
observe that since r" is bounded, 

<C\D^u]\', 



and similarly {D*^RJ)^ < |_D*^w"| . Defining 
Ax 



By the assumption (4.10) G" and q2 are locally bounded, hence we have that 
\elj\<CAx\D^u]\'^, fori = 3, 4, 

and for some C which is independent of Ax. By the same argument used to prove 
Lemma 4.2, also e^^Ax and e^^Ax are in Aiiod^)- We now have established that 

D^it^+D^q^^=EZ^, 1 = 1,2, 

where 

ill^ = - k) , rjl^ = G{R]) - G{k), 
and q^j are given by the corresponding fluxes 

ql^ = G{Ri.i) - G{k), qi^ = G'{ef de. 

Jk 

The terms E"j are such that the corresponding piecewise continuous functions 
Ei,Ax are in M\oc{^)- 
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Lemma 4.3. Suppose (4.10) and the assumptions A.l - A. 4 hold. Then there 
is a subsequence of {-Rax} (i^-ot relabeled) and a function R such that Rax ^ R 
a.e. {x,t) e n. We have that R £ L°°{[0,T]; L^{R)). Furthermore, R satisfies 

jRt + G{R).^ < 0, xeR, t>0, 
\R= |^io|^ xeR, t = 0, 

in the distributional sense. 

Proof. The condition A. 3 means that if we can show that (rii^Ax)t + {<li,Ax)x for 
i = 1, 2 is compact in H^^{^^), Rax converges strongly. To show this, let V' be a 
smooth function in Hl^^{il). We have that 

{{'ni,Ax)t + {qiAx)x, V') = X! / / ^'"i'^* + Ci^^ ^^^^ 

= / ^"^ - ijj (t„, a;)) dx 

/tn+l 
llj {i' {X] + l/2,t) - i) {Xj-.i/2,t)) dt 

raj •'Xj^l/l 

= -y2{Dlv?j+D-q^j) II ^{x,t)dxdt~ I T^,^Ax{x,0)tp{x,0)dx 

+ E J I ^-(2;, - (tn+i, a;) dxdt 

n,j ^7 

= - ^ E^j J J ijj{x, t) dxdt - j rj.^Axix, O)'0(a;, 0) dx + ai + a2- 

n.j * 

Hence, {r]i^Ax)t + {qi,Ax)x consists of the sum of ai and a2 and a term which is in 
■M\oc{^)- We must show that ai and a2 are compact in iJ^^(ri). First, observe 
that since i?" is uniformly bounded, 

{D'+vIjY < C and [D^ql^f < C \D^u]\'' . 

Using this 

ai < \D-u]\ JJ J"^' \Mx,t)\ dr dxdt 

n,j ^7 * 

<Cy^\D^u^\ II ^t,,+i-t( I "^\Mx,T)f drV^^dxdt 
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1/2 



dxdt I 



<C^t(AtAx^Ax\D_vJ^f\ ' (S^ 11 {iPt{x,t)f dxdt\ ' 
(4.11) <C%/A^||V||^i(o)- 

Hence ai is compact in H^^{ fl), that 02 is compact follows by analogous arguments. 
Now we have established the strong convergence (along a subsequence which we do 
not relabel) of i?Ax- 

By the CFL-condition, (4.7), the right hand side of (4.8) is negative, hence 

Multiplying this with a non-negative test function ip, doing a summation by parts 
and then sending Ax to zero, using the same arguments that led to (4.11), yields 
that the limit i? is a subsolution in the distributional sense. □ 



Since Rax converges strongly to R, also r^x will converge strongly to r :— \/R. 
The sequence {uax} a^^q is uniformly bounded, so a subsequence will converge 
weak-* to some function u £ L°°{il). By using the arguments leading up to (4.11) 
it is straightforward to show that 

wpt + u(f>{r)^px dxdt + / uq{x)^P{x,0) dx — 0, 

for all i/j € C^(ri). Hence the limit u is a distributional solution of 

ut + {u(l){r))^ = 0. 

In order to conclude that u is a weak solution to (1-1), we would have to show that 
|u| = r. We have not been able to prove this, and merely conclude that |u| < r. 
The reason for this is that v i-> \v\ is convex, and that weak limits of a convex 
function are not less than the convex function of the weak limit. 

To overcome this difficulty, we propose another fully discrete scheme based on 
explicit decoupling of the variables r and w. 

4.1. A scheme which enforces the entropy condition. Define 

if VAx 7^ 0, 
if rAx = 0, 



AtD^ff, n>0. 




At(l)]D^w^, n>0, 



To ensure the convergence of the approximations {r^x} we choose 

(4.14) Ai||/'||^^(K) < Ax. 

We list some useful properties of r^x in the next lemma [■)]. 
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Lemma 4.4. Assume that the CFL condition (4.14) holds and ro € BV{M.) n 
i°°(M). Then for each Ax > we have that 

(a) —M < r/^xix,t) < M, for all x and t > 0. 
(h ) For n > the functions 

are non-increasing. In particular this means that the family {r'^^}^ 
(uniformly in Ax) hounded m L°° {M.+ ; {R)) n BV(M x M+). 
(c) Moreover r^xi'it) — ^ i^i'^t) strongly in L^{M.) for all t > 0, where r G 
Lip([0,r];Li(M)) and is the unique entropy (in the sense of Kruzkov) so- 
lution of the conservation law 



(4.15) 



|rt + /(r), =0, 
|r(a;,0) = tq. 

Observe that, we can write the scheme (4.13) as 



= (1 - \^^) + A0>;_i. 
If A0" < 1 for all j, then wj~^^ is a convex combination of w" and w"„]^. Thus 

(4.16) infw° < u;^" < supw°, 7i > 0. 

^ j 

and 

(4.17) |u>Aa;(-,t)lB.v.(R),. < \wAxi-,0)\B,v.{R)^ ■ 

Furthermore 
(4.18) 

AxY,\w]+' -Wj\< At\m^^ K - < CAt\wAA;tn)\B.V.iM)^ ■ 

j 3 

Hence the map t i— WAxi'^t) is L^- Lipschitz continuous. Finally, the above es- 
timates (4.17), (4.18) and an application of Kolmogorov's compactness criterion 
(Lemma 2.3) shows that w = limA^^o "Waxj where the convergence is along a sub- 
sequence. Furthermore w e C([0,r]; (^^^^(M))"). 

AtD^f-) {w--At4>-D.u,-) 



Multiply the equation (4.12) for r"+^ with that (4.13) for -u;"+^ to get 



r>" - At (w'W.f'' + f''D_w'l) At^^'lD^f^'D^w 



>" - At + - At {f^ f^^ D_w- 



+ At'^(f>]D^f^D^w^ 
= r>^" - AtD^ (f^wj) + At (/; - D^w^ {X^] - 

Then we have 

(4.19) Dl (r>;.0 -f D_ (/>;.0 = At [Xc^] - l) (/; - D_w- =: e,' , 

Let now tp G C^(f2) be a test function, multiply the above equation by ^ and 
integrate over O to get 

Y / r^w^Dt'ip + f;i'w^D+tPdxdt 
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Since we have the convergence of r^^ and wax , the left hand side of this converges 
to 



rwipt + fif)ui4'x dxdt + j r{x,0)w(x,0)'il}{x,0) dx. 
Regarding the right hand side we have 

^ / / e^ipdxdt 



n,j 



< At Uh^^a) (A + 1) At 5] |/; - /j^il 



< AtCllV'llioo(K) ||wAx(a;,0)||^„o(jj)„T||rA:E(a;,0)||s^^(jj) , 
where T is such that suppV' C [0,T]. Hence 

rwipt + fi''^)w'4'x dxdt + / r{x,0)w{x,0)^{x,0) dx — 0. 

Hence, we see that rw is a weak solution to the Cauchy problem 

irw)t + {f{r)w)^ = 0. 

In other words, (r, rw) is a weak solution to 

rt+f(r)x = 0, 
(4.20) t J \ jx 

{rw)t + {cf){r)rw)x = 0. 

Next, we shall make use of Lemma 2.4 to conclude that \'w\ = 1. 

Finally, collecting all the results above, we have proved the following theorem. 

Theorem 4.1. Assume that uq e B.V.{R). If X ^ At/ Ax satisfies the CFL- 
condition A < sup^. /'(|uo(x)|), and uax is defined by (4.12), (4.13), then u = 
liniAs-i-o '^Aa; the unique entropy (in the sense of Definition 1.2) solution to (1-1). 

Remark 4.1. Observe that to establish the convergence of the entire sequence 
{u£^^x\ Ax>Q' /i?^st use Kolmogorov's compactness criterion to establish conver- 
gence of a subsequence, then, since the entropy solution u is unique, every subse- 
quence o/ {wAa:} will Contain a subsequence converging to u. This means that the 
entire sequence converges. 

We propose another scheme based on discretizing the "conservative" form (1.6)- 
(1.7). Let TAz: and max satisfy 



(4.21) 



r-0 — U/0| 



AtD_f{r^), n>0. 



and 

(4.22) u^+' = ~ AtD^ {u](j,{rj)) , 

for n > and /(r) = r(j){r), with u° given by (3.2). Again as before, regarding the 
convergence of the approximations {rAx} we have Lemma 4.4. On the other hand, 
we can rewrite the scheme for as 
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Each component = F" 



\ ) is linear and monotone in both 



its arguments. Furthermore, the scheme for r"^ can be written 



Set 



Then, by the hnearity of FJ", 



^-1 



Therefore inf ., w° < w"^^ < sup,- w?. We have that 

■> J J J J J 



'uf'''^ > 0, therefore 



(4.23) 



\uAx{x,t)\ < rAx{x,t) < \uo{x)\ . 



Hence, the sequence {uax}ax>o uniformly bounded, and there converges 
weak-* (modulo a subsequence) to some function u. A straightforward computation 
yields that 



UAx'4't+'UAx(f'{rAx)4'tdxdt+ I Ua{x)'ijj{x,0)dx 

n 

< E 

n=l,j 



I? 



_ ^j) dxdt +\u]\(l) (r^) / / {D+i; - ^j^) dxdt 



1 



At 



UAx{x, Q)i!{x, t) — lS.tuQ{x)%l)(x, 0) dx 



< CAx f ||V'tt||i==(o) + W^^xxWl^^u) + ll'^tllL~(n) 



Sending Ax to zero, we see that the limit u is a weak solution to the Cauchy 
problem 

ut + {u(l){r))^ = 0. 

Again we shall make use of Lemma 2.4 to conclude that |u| = r. 
Finally, We have proved the following theorem. 

Theorem 4.2. Assume that uq e B.y.(M). If X ^ At/Ax satisfies the CFL- 
condition A < sup^, /'(|iio(x)|), and uax is defined by (4.21), (4.22), then u = 
liniAs-i-o '^Aa; the unique entropy (in the sense of Definition 1.2) solution to (1-1). 



5. Numerical experiments 

We close this paper by demonstrating how these schemes work in practice. We 
perform all the computations for 2x2 system with 0(r) = r^. 

5.1. Numerical Experiment 1. In this case we approximate the system (1.1) 
with initial data 



(5.1) 



C/o(x) 



a: < 0, 
X > 0. 



It is not difficult to find the exact solution of (1.1) in this case. For the sake of 
completeness we write the expHcit form of the exact solutions U{x,t) — tj{x/t). 
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If \Ui\ < \Ur\, then 
(5.2) [7(0 
If \Ui\ > \Ur\, then 



m'<^<3\Ui\\ 
(|)V2|^, 3|C/,|^<e<3|C/.|^ 



2 ^ /■ ^ iT-r |2 , ITT I ITT \ , Irr |2 



UiO = { Urn, |f//r < e < l?:^;! + \Ul\ \Ur\ + \U, 

^>\Ulf + \Ul\\Ur\ + \Ur\\ 

with Urn — ^f\Ur in both cases. 

In what foUows, we test the fuUy discrete exphcit numerical scheme (4.1) with 
initial data 

fc/_, x<0, 
X > 0, 

where 

[/_ = (0.5, 1.5), C/+ = (1.5,2.0), 

for the first experiment and 

= (1.5,2.0), [/+ = (0.5,1.5), 

for the second experiment. The computations are performed on a computational 
domain [—1,20] with 4000 mesh points. To enforce the CFL condition we set the 
time step At = (C F L) Ax / 3 sup \Uo\ , where we use a CFL number 0.75. Although 
we do not plot the exact solutions, a comparison of the computational results dis- 
played in Figs 5.1 with the exact solution shows good agreement. 

5.2. Numerical experiment 2. In this case, we test our fully discrete explicit 
numerical scheme (4.12)-(4.13) with initial data Uq — tq-wq, where 



with 



a; < 0, 
a; > 0, 



r_ = 1.0, r+ ^ 0.75, 
for the first and third numerical experiments and 

r_ = 0.75, r+ = 1.0, 
for the second and fourth numerical experiments. Similarly, for wq we take 

{(1.0,0.0), a;<0.2 

(cos(87r(x - 0.2)), sin(87r(a; - 0.2))), 0.2 < x < 0.7, 

(1.0,0.0), a;>0.7, 
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(a) T = 



(b) T = 




(c) T = 0.50 (d) T = 0.50 

Figure 5.1. Left column: Experiment- 1: The dotted-dashed 
curve represents the first component of U, the solid curve repre- 
sents the second component. Right column: Experiment-2: The 
dotted-dashed curve represents the first component of U, the solid 
curve represents the second component. 



for the first and second numerical experiments and 



Wo[X) = 



(1.0,0.0), X < 0.2, 

(-1.0,0.0), x>0.2. 



In this case also, it is easy to find the exact solution. Although we do not plot the 
exact solutions, we give the explicit form of the exact solution. The exact solution 
is given hy U — rw with 

'■[X, t) = <. with s = r_ + r_r+ + r,, 

X > st, 
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and 

wo{^{x — r'^t)), r'^t<x<st, 
Wq{x — r^t), X > st, 

for the first and third numerical experiments. Similarly, 

x<3rH, 

r{x,t) = {{x/3t)^/^, 3rH<x<3rlt 



T+, X > 3r^i, 

Wo{x — r'^t), X < r^t, 

WaC^ix -rlt)), rlt <x <'ir'it, 



and 



w{x, t) 



Wq{x — r'^t), X > 3r^<, 

for the second and fourth numerical experiments. 

In all the experiments computational domain is [—1,4] and we use Neumann 
boundary conditions at the left boundary. We also use a CFL number 0.75 and 
4000 mesh points for all the experiments. A comparison of the computational 
results displayed in Figs 5.2-5.3 with the exact solution shows good agreement. 

Below we show the computational results for four different qualitatively signifi- 
cant sets of data: a compression or an expansion wave in r initiated slightly behind a 
continuous pulse or a discontinuous contact wave in w. Fig 5.2-5.3 display the com- 
puted solution at three different times. In the plots, the dot-dash curve represents 
the first component of U and the dotted curve represents the second component, 
while the solid curve represents the r-component of (r, U). 

5.3. Numerical convergence rates. We have not obtained any theoretical con- 
vergence rates for the schemes presented here. Never the less, it is interesting to 
check the possible convergence rate in practice. To this end we have used Riemann 
initial data (5.1) with Ui = (1, 1) and C/,- = (3, 1). In this case the exact solution is 
given by formula (5.2), i.e., a rarefaction wave followed by a contact discontinuity. 
We define the relative error for a scheme as 



(5.3) £; = 100 X 



j:^\uf-u{x„NAt)\ 
j:Au{x„NAt)\ 



where u is the exact solution found by (5.2) and is the approximation computed 
by the numerical scheme. Note that this is a first order accurate approximation to 
the relative error since u is piecewise continuous. 

We have computed the errors for the three schemes (4.1) ( "schemel" ), the conser- 
vative scheme (4.21) - (4.22) ("scheme2") and the non- conservative scheme (4.12) 
- (4.13) ("schemes"). Table 5.1 summarizes the results. We computed the ap- 
proximations for < = 1, and used Ax = 40/2^ for = 5, . . . , 14 for x e [—1, 39]. 
From this table it emerges that the three schemes produce very similar errors and 
numerical convergence rates. This convergence rate is expected to be not higher 
than 1/2, since the solution contains a contact discontinuity, and the schemes are 
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(a) T = 



(b) T = 




(c) T = 0.25 



(d) T = 0.25 



(e) T = 0.75 



(f) T = 0.75 



Figure 5.2. Left column: Experiment- 1: A shock wave initiated 
behind a continuous rotational wave. The dotted-dashed curve rep- 
resents the first component of C/, the dotted curve represents the 
second component and the solid curve represents r. Right column: 
Experiment-2: An expansion wave initiated behind a continuous 
rotational wave. The dotted-dashed curve represents the first com- 
ponent of U, the dotted curve represents the second component and 
the solid curve represents r. 
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(a) T = 



(b) T = 



(c) T = 0.25 





/ 




k 


(d) T = 0.25 




(e) T = 0.75 



(f) T = 0.75 



Figure 5.3. Left column: Experiment-S: A shock wave initi- 
ated behind a discontinuous rotational wave. The dotted-dashed 
curve represents the first component of J7, the dotted curve repre- 
sents the second component and the solid curve represents r. Right 
column: Experiment-4: An expansion wave initiated behind a dis- 
continuous rotational wave. The dotted-dashed curve represents 
the first component of U , the dotted curve represents the second 
component and the solid curve represents r. 
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scheme 1 
E rate 



scheme2 
E rate 



schemes 
E rate 



N_ 
"5" 
6 
7 
8 
9 
10 
11 
12 
13 
14 



3.32 
2.04 
1.31 
0.81 
0.51 
0.32 
0.20 
0.13 
0.09 
0.06 



0.70 
0.64 
0.69 
0.67 
0.68 
0.64 
0.65 
0.62 
0.60 



3.36 
2.08 
1.35 
0.83 
0.52 
0.32 
0.20 
0.13 
0.08 
0.05 



0.68 
0.63 
0.69 
0.69 
0.68 
0.67 
0.67 
0.65 
0.63 



3.40 
2.31 
1.50 
0.89 
0.54 
0.33 
0.21 
0.13 
0.08 
0.05 



0.55 
0.62 
0.75 
0.72 
0.70 
0.68 
0.68 
0.65 
0.63 



Table 5.1. Relative errors and rates. 



formally first order. Table 5.1 actually seems to predict that the three schemes 
have a convergence rate of about 0.6. 
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